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ABSTRACT 

Gamma ray burst afterglows can be identified in single epoch observations using three or more 
optical filters. This method relies on color measurements to distinguish the power law spectrum 
of an afterglow from the curved spectra of stars. Observations in a fourth filter will further 
distinguish between afterglows and most galaxies up to redshifts z « 1. Many afterglows can also 
be identified with fewer filters using ultraviolet excess, infrared excess, or Lyman break techniques. 
By allowing faster identification of gamma ray burst afterglows, these color methods will increase 
the fraction of bursts for which optical spectroscopy and other narrow-field observations can be 
obtained. Because quasar colors can match those of afterglows, the maximum error box size where 
an unambiguous identification can be expected is set by the flux limit of the afterglow search and 
the quasar number-flux relation. For currently typical error boxes (10 - 100 square arcminutes), 
little contamination is expected at magnitudes R < 21.5 ± 0.5. Archival data demonstrates that 
the afterglow of GRB 000301C could have been identified using this method. In addition to 
finding gamma ray burst counterparts, this method will have applications in "orphan afterglow" 
searches used to constrain gamma ray burst coUimation. 



Subject headings: Gamma rays — bursts 
1. Introduction 

The search for counterparts to gamma ray 
bursts (GRBs) at longer wavelengths has been 
recognized from the beginning as a key to un- 
derstanding the bursts' nature. This motivated 
models that predicted gamma ray burst after- 
glows from hot material in relativistically expand- 
ing GRB remnants (Paczyhski & Rhoads 1993; 
Katz 1994; Meszaros & Rees 1997). The afterglow 
was predicted to arise from synchrotron radiation, 
and its spectrum is therefore well described by a 
broken power law. Both the prediction of fading 
counterparts with power law spectra and the gen- 
eral expectation that GRB counterparts would tell 
us a great deal about the bursts have been amply 
borne out by observation since the first successful 
observations of GRB afterglows in 1997 (Costa et 
al 1997; van Paradijs et al 1997; Frail et al 1997). 

The greatest rewards so far have come from 
detailed studies of individual afterglows with the 
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best time and spectral coverage. Rapid after- 
glow identifications are a key to achieving such 
data sets, because some essential types of data 
cannot be obtained before the afterglow position 
is known to 1 arcminute (near-infrared and 
submillimeter photometry) or even ~ 1 arcsec- 
ond (spectroscopy). Traditionally, afterglows are 
found through a variability search within GRB er- 
ror boxes. This requires multiple epoch obser- 
vations whose separation in time At must be a 
reasonable fraction (> 0.3) of the time elapsed 
from the GRB to the first epoch. This lag is usu- 
ally further increased by telescope scheduling is- 
sues to ~ 1 day. Moreover, afterglow light curves 
occasionally include "plateau" phases (e.g. GRB 
970508 ait <2 days [Pedersen et al 1998]; GRB 
971214 at t < 0.3 days [Gorosabel et al 1998]; and 
GRB 000301C at t ~ 6 days [Bernabei et al 2000; 
Masetti et al 2000]). Two-epoch observations of 
such afterglows could miss them entirely. 

This delay in afterglow identification is not nec- 
essary, because afterglows are characterized not 
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only by fading behavior but also by power law 
spectra. Such spectra are relatively rare in the 
optical sky, which is dominated by Galactic stars 
at bright flux levels and by galaxies at faint flux 
levels. Relative to a stellar spectrum, a GRB af- 
terglow is blue at blue wavelengths (< 0.4/zrn), 
and red at red wavelengths (> O.G/im). A color- 
color plot should therefore isolate a GRB afterglow 
from nearby stars efficiently, provided it spans a 
sufficient wavelength range to detect spectral cur- 
vature (or the Balmer jump) in stars of any sur- 
face temperature. Under many circumstances it 
will also be possible to identify afterglows using 
simpler though less robust single color criteria (ei- 
ther ultraviolet excess ["UVX"] or infrared excess 
["IRX"]). Moreover, absolute photometric calibra- 
tion is not required so long as a sufficient solid an- 
gle is observed to empirically determine the stellar 
locus in color space. 

Section 2 describes model color-color plots that 
show the expected loci of GRB afterglows, stars, 
galaxies, and quasars. Section 3 demonstrates em- 
pirically that the color-color method could have 
identified the afterglow of GRB 000301C in a sin- 
gle epoch. Finally, section 4 considers advantages 
and limitations of the method, including confu- 
sion limits as a function of flux level, likely se- 
lection effects in the sample of afterglows found 
through color-based searches, and applicability of 
the method to "orphan afterglow" searches. 

2. Color-Color Space Models 

To explore this method in detail, I have calcu- 
lated synthetic color space locations for GRB af- 
terglows, stars, galaxies, and quasars. Model cal- 
culations were performed for the widely available 
Johnson U, B, and V filters and Cousins R and 
I filters, though color-color searches will work for 
any set of filters spanning a sufficient wavelength 
range. Filter transmission curves were taken from 
Bessell (1990). The photometric zero point for all 
Johnson-Cousins system synthetic colors was set 
to the Vega model atmosphere of Kurucz (1979). 

Model spectra of objects at high redshift were 
attenuated using the mean intergalactic absorp- 
tion given by Madau (1995)^. For redshifts z>2, 

^The photoelectric absorption was treated using the approx- 
imation to Madau's equation 16 from his footnote 3, which 
is within 5% of the full calculation. 



the flux in the bluest fllter can be greatly reduced 
by neutral hydrogen absorption in the intergalac- 
tic medium. This absorption consists of a super- 
position of Lyman a forest lines, higher order Ly- 
man lines, and bound-free Lyman continuum ab- 
sorption (Madau 1995). The net effect is a shift in 
the bluest color measured. The relevant redshift 
depends on the choice of bluest filter. The Lyman 
a line is redshifted to the peak transmission wave- 
length of the Johnson U band at z = 2.04, and the 
Johnson B filter around z = 2.5. 

In addition, I calculated reddening vectors for 
dust both in the Milky Way and in afterglow host 
galaxies at a range of redshifts, using the empiri- 
cal fits of Pei (1992) to the extinction laws of the 
Milky Way and the Magellanic Clouds. 

In the following subsections, I describe the color 
calculations for each class of object considered, 
and then summarize the results. 

Stars: Stellar colors were derived from the 
spectrophotometric catalog of Gunn and Stryker 
(1983), which spans spectral types from O to M 
and luminosity classes from the main sequence to 
supergiants. 

GRB Afterglows: Afterglows were taken as 
pure power law spectra, f^, oc j^", with indices 
0.5 < —a < 1.5. Such spectra describe syn- 
chrotron emission from a power law distribution 
of electrons away from break frequencies. Devia- 
tions from power law spectra occur near spectral 
breaks (see Sari, Piran, & Narayan 1998) and in 
the case of interstellar absorption (section 4). The 
range of spectral indices considered corresponds to 
the commonly measured range of electron energy 
indices, 2 < p < 3. Here p describes the number- 
energy distribution of electrons accelerated at the 
GRB remnant's external shock, N{£) oc £~p, and 
I have considered spectra both above and below 
the cooling frequency. 

Galaxies: Galaxy spectra were taken from 
the GISSEL stellar population synthesis code 
(Bruzual & Chariot 1993) for a few representative 
star formation histories. No attempt was made to 
include galaxy evolution effects, wliicli may result 
in elliptical galaxy colors that are slightly too red 
at high redshift (see, e.g., Bruzual [1996]). The 
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resulting galaxy colors agree with those computed 
by Fukugita, Shimasaku, & Ichikawa (1995) to 
within 0.1 mag for most filters, galaxy spectral 
types, and redshifts; this is acceptable given the 
differences in the input spectra in the two works. 

Galaxy spectra at optical wavelengths are a 
superposition of diverse stellar spectra, modi- 
fied slightly by nebular emission and dust effects. 
This superposition implies that galaxies have less 
sharply peaked spectra than individual stars and 
so lie on the same side of the color-color space stel- 
lar locus as power law spectra. Indeed, the spec- 
trum of a star forming galaxy can be a good power 
law over most of the optical window. Fortunately, 
the 4000A break due to stellar atmospheric fea- 
tures interrupts this approximate power law and 
can be used to distinguish almost all galaxies at 
z < 1 from afterglow candidates given sufficient 
color information. With observations in three fil- 
ters, there is some redshift where the 4000A break 
lies in the middle filter, and galaxy colors can 
match those of a pure power law. The addition 
of a fourth filter breaks this degeneracy. For ex- 
ample, in the U — B, R — I plane, an irregular 
galaxy at z « 0.4 lies on the power law sequence; 
however, its B — R color is much redder than the 
corresponding power law. In addition to color cri- 
teria, it is of course possible to distinguish galaxies 
from GRB afterglows morphologically. The frac- 
tion of galaxies that appear extended depends on 
the available spatial resolution, with the a few 
galaxies appearing pointlike even in Hubble Space 
Telescope images, but subarcsecond resolution will 
in practice greatly reduce confusion due to galax- 
ies. 

Quasars: Quasar colors were calculated using 

the composite qiiasar spectrum by Francis et al 
(1991), extrapolated to A > 0.6/xm with fx oc 
^-2.18 rjij^-g extrapolation will result in modest 
color errors for low redshifts where the Ha line 
falls in a measured filter, but should not affect our 
conclusions greatly. 

Quasar spectra are reasonably well approxi- 
mated by power laws (plus broad emission lines) 
across a wide wavelength range. Color-based cri- 
teria have been used to identify quasars for nearly 
40 years (see the review by Burbidge 1967), and 
the use of far-red filters to supplement bluer UBV 
passbands was suggested as early as 1967 (Brac- 



ccsi 1967; Braccesi, Lynds, & Sandage 1968). 
Quasars are therefore extremely difficult to dis- 
tinguish from GRB afterglows using only broad 
band colors. The main discriminant available is 
that quasars tend to be somewhat bluer than typ- 
ical afterglows. In practice, we expect confusion 
with quasars to set the maximum area over which 
an unambiguous afterglow identification can be ex- 
pected. This confusion can be resolved with low 
resohition spectra, which identify most quasars 
unambiguously by their strong, broad emission 
lines. 

2.1. Results 

Sample color-color diagrams are shown for the 
U — V, V — I plane in figure 1, and the B — V, R—I 
plane in figure 2. These have been chosen as repre- 
sentative cases where the two-color method works. 
Additional combinations of filters are also effective 
provided that at least three filters are used, the 
bluest filter is at least as blue as the Johnson B 
band, and the reddest is at least as red as Cousins 
I. Unfortunately, the observationally easiest filters 
(BVR) do not yield an adequate separation be- 
tween afterglow and stellar colors to expect reli- 
able identification of afterglows, barring unusually 
high photometric precision. 

Given a sufficient filter set, the afterglows are 
well separated from the stellar locus. Moreover, 
the reddening vector for Milky Way dust runs 
essentially parallel to both the stellar sequence 
and the power law spectral sequence. Thus, even 
though reddened afterglow spectra are no longer 
strict power laws, they remain distinct in color- 
color space. 

Of greater concern is dust in the GRB host 
galaxy. Because this dust is at high redshift, 
its 2 175 A feature can enter the observed wave- 
length range, changing the direction of the red- 
dening vector in color-color space. For redshifts 
2.2 ^ z < 3.3, this feature falls in the reddest 
filter of the color-color plot, thus making the long- 
wavelength color (e.g., V-I) slightly bluer. A mod- 
erate amount of dust (Ay = 0.5) can then move 
the afterglow locus onto the stellar locus in color- 
color space. However, this effect is only impor- 
tant for Milky- Way type dust. The most impor- 
tant consequence of dust for color-based afterglow 
searches is therefore not reddening but extinction, 
which can reduce the received flux below the de- 
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tcction limit of the search whatever the reddening 
law. However, GRBs can destroy dust at distances 
up to 10 pc (Waxman & Draine 2000) or beyond 
(Fruchter, Krolik, & Rhoads 2001), reducing these 
concerns substantially for bursts at high Galactic 
latitude. 
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Fig. 1. — Model U-V, V-I color-color plane. 5- 
pointed stars mark the location of GRB after- 
glows, represented by pure power law spectra with 
indices fi, oc i^" for a = —0.5, —1.0, —1.5. Small 
points represent stars from the Gunn-Stryker 
(1983) atlas. Dotted lines show the locations of 
different galaxy spectra as a function of redshift. 
Star formation activity increases from lower right 
(Elliptical galaxy model) through spiral models 
(Sb and Sc) to an actively star forming Magel- 
lanic irregular (Im) model. Open circles along 
these tracks mark intervals of 0.5 in redshift, from 
zero to 3.0. No galaxy evolution is incorporated in 
these tracks. Finally, the solid line shows the locus 
of quasars as a function of redshift. Filled points 
along this track again mark redshift intervals of 
0.5 from zero to 3.0. Intergalactic absorption will 
shift the afterglow locus to redder U — V at fixed 
y - / for z > 2.0. 



3. Empirical Tests 

Archival multicolor data is available for a few 
GRB fields from the U.S. Naval Observatory's 
GRB foUowup team (see Henden et al 2000). I 
have used this data for the field of GRB 000301C, 
together with published photometry of the after- 
glow (Jensen et al 2000), for an empirical demon- 
stration of color-based optical afterglow searches. 

Two color-color plots for this field are shown 
in figure 3. In both cases, we see that the af- 
terglow is an outlier in color-color space. In the 
U — V, V — I figure, there are a few other outliers, 
all of which are systematically redder than the af- 
terglow in both colors. In the B — V, R — I plot, 
the afterglow is the single most dramatic outlier 
and would be the first target chosen for prompt 
spectroscopy or other rapid narrow-field followup 
observations. 

4. Discussion 

Color-based searches have the potential to iden- 
tify GRB afterglows faster than a two-epoch vari- 
ability search, and hence to allow more thorough 
followup of many afterglows. There are three ba- 
sic color selection criteria available for afterglows: 
IR excess ("IRX"), UV excess ("UVX"), and Ly- 
man break. The appropriate combination of these 
depends on both the redshift of the burst and 
the available instrumentation. The primary fo- 
cus of this paper is on color-color plane methods, 
which combine the IR excess and UV excess cri- 
teria and are therefore more robust than either 
criterion alone. 

Still, as the observations of GRB 000301C 
demonstrate, potentially occupied regions of color- 
color space remain vacant in practice for high- 
latitude fields observed to a limiting magnitude 
V ^ 21. In particular, there are no very hot stars, 
few or no high redshift galaxies, and few or no 
quasars. Under these circumstances, a simple UV 
Excess ("UVX") test could identify the afterglow 
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Fig. 2. — Model B-V, R-I color-color plane. Sym- 
bols are the same as for figure 1. Because neutral 
hydrogen absorption first becomes important at 
higher redshift in the B band than the U band, 
color tracks for galaxies and quasars extend to 
somewhat higher redshift here, and the afterglow 
locus shifts to redder B — V for z > 2.5. 




Fig. 3.— Color-color plots for the field of GRB 
000301C. Data are taken from Jensen et al (2000) 
for the GRB afterglow (marked by a five-pointed 
star symbol); and from Henden et al (2000) for 
other objects in the field. For the Henden et al 
data, error bars on U-V were computed by adding 
in quadrature the reported U-B and B-V errors; 
similarly, V-I errors were computed from V-R and 
R-I. The GRB is an outlier in both plots, and the 
most dramatic outlier in the B-V, R-I plane. 



on the basis of a single color measurement (prefer- 
ably including one filter at A < 4000 A). 

Likewise, at near-infrared wavelengths, most 
stars and galaxies will have spectral slopes near 
the Rayleigh-Jeans value (2) and will therefore 
be bluer than afterglows. This is again demon- 
strated by GRB 000301C, for which we found 
J ~ K' ^ 1.45 (Rhoads & Fruchter 2001), about 
0.3 magnitudes redder than even the coolest stars 
tabulated by Johnson (1966). However, UVX and 
IRX methods are vulnerable to the presence of 
very hot stars and highly reddened sources respec- 
tively, while the two color method is not. 

The Lyman break (or "dropout") method uses 
intergalactic absorption to find high redshift after- 
glows whose strongly attenuated blue flux places 
them in a region of color space occupied by es- 
sentially all classes of high redshift object, on the 
opposite side of the stellar locus from low red- 
shift afterglows. The observational requirements 
to find these bursts from Lyman break colors are 
more exacting. Long integrations in the bluest fil- 
ter are required to conclusively identify "dropout" 
afterglows, while lower redshift afterglows are un- 
usually bright in the bluest filter and require cor- 
respondingly shorter integrations. Even so, most 
afterglows are brighter than their host galaxies at 
discovery, which means that a high redshift after- 
glow may be an easier dropout to identify than a 
high redshift galaxy. 

Color-based searches are subject to certain se- 
lection effects, which must be considered in inter- 
preting their results. These effects include: (1) 
Redshift: The onset of intergalactic hydrogen ab- 
sorption in the bluest filter observed determines 
the highest redshift where UV excess methods can 
work and the lowest redshift where Lyman break 
methods can work. The UVX method will fail at 
z « 2 for U band and z w 2.5 for B band, while 
the Lyman break method first becomes practical 
in the U band at z « 3.3. Thus, two-color opti- 
cal searches will not identify afterglows between 
about z — 2.5 and z — 3.3. The IR excess will 
continue to be practical to quite high redshift, e.g. 
to 2 ~ 15 for a J — iiT color, since the onset of neu- 
tral hydrogen absorption simply accentuates an IR 
excess until the absorption reaches the redder fil- 
ter. Because there is considerable variance in the 
hydrogen column density along different lines of 
sight (Madau 1995), the redshift boimdaries dis- 
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cussed here are not sharp demarcations. Rather, 
the detection efficiency declines smoothly, and the 
characteristic redshifts given above correspond to 
about 50% of peak efficiency. 

(2) Heavily obscured afterglows can drop be- 
low the search's detection threshold. Extinction 
affects any search method, though it is exacer- 
bated for UVX and two-color methods by the 
increase of extinction towards bluer optical-UV 
wavelengths and the nc;c;cssity of observations in 
observer frame blue light. Additionally, there will 
be a mild selection against finding afterglows in 
galaxies with Milky Way type dust extinction at 
redshifts 2.2 < z < 3.3 using two-color optical 
searches. However, this overlaps the stronger se- 
lection due to intergalactic Hydrogen absorption 
to such an extent as to be almost unimportant. 

(3) For gamma ray bursts occurring near dense 
molecular gas, absorption by excited molecular hy- 
drogen can occur for rest wavelengths as long as 
1650A (Draine 2000). Absorption between 1650 
and 1300A can reach 50% for H2 column densities 
around 10^^ cm^^. This will act to exclude after- 
glows in molecular clouds from UVX or two-color 
samples at redshifts beyond 1.3 (using U band) or 
1.65 (using B). 

(4) Because afterglow spectra are broken power 
laws rather than pure power laws, there will be 
times in the evolution of each afterglow when its 
colors deviate from the models in section 2. Breaks 
in synchrotron spectra are not absolutely sharp, 
so this is unlikely to be a major selection effect 
for purely optical searches, where the wavelength 
range spanned is only a factor of ~ 2.5. Searches 
employing larger wavelength coverage (UV or IR) 
will be more sensitive to this effect. 

(5) Afterglows located atop brighter host galax- 
ies may be hard to find because the dominant color 
is that of the galaxy. This effect will depend on 
instrumental spatial resolution as well as the in- 
trinsic properties of the afterglow and host galaxy. 
(Variability searches can be similarly compromised 
unless image subtraction methods are used.) 

(6) An advantage of this method is that short- 
term "plateaus" in afterglow light curves do not 
affect it. It may therefore find some afterglows 
that traditional variability searches miss. 

Even afterglows that lie in the predicted region 
of color space may be difficult to identify by these 



methods in cases where the error box is so large 
that a substantial number of other sources occupy 
the same region of color space. The critical error 
box size where one "confusing" source is expected 
depends on the magnitude limit of the search. Ul- 
timately, the quasar number-flux relation gives a 
lower limit to the confusion level. Additionally, 
for three filter searches, star forming galaxies at 
^; ~ 0.5 can add to the confusion. To achieve a 
unique afterglow identification, it is necessary to 
match the magnitude limit of the search to the size 
of the error box, and to observe the field while 
the afterglow remains brighter than this magni- 
tude limit. A unique identification may not al- 
ways be needed, however. In particular, a list of 
several viable candidates may be enough to obtain 
a same-night afterglow spectrum using a multi- 
object spectrograph with a sufficient field of view 
and rapid setup procedure. The correct candidate 
could then be identified later based on either spec- 
tra or variability. 

The quasar number-magnitude relation for z < 
2.2 and magnitude B < 24 is predicted by Malho- 
tra & Turner (1995) for both {n = 1, A = 0) and 
(fi = 0.1, A — 0.9) cosmologies, based on observa- 
tions to S < 22 by Boyle et al (1987) and Boyle, 
Shanks, & Peterson (1988). 

The approximate range of redshift where star 
forming galaxies overlap the GRB afterglow locus 
in color-color space is 0.35 ^ z < 0.6 for the U — V, 
V - I diagram, or 0.45 < -2 < 0.6 for the B -V, 
R — I diagram. I estimate a number-magnitude 
relation for such galaxies using both the local lu- 
minosity function from Loveday et al (1992) and 
direct covmts of galaxies with 0.55 < y — / < 1.28 
and 0.4 < z < 0.6 from the Canada-France Red- 
shift Survey (Lilly et al 1995a; Le Fevre et al 1995; 
Hammer et al 1995). The predicted number counts 
amount to approximately 10% of all galaxy rnim- 
bcr counts (see, e.g., Lilly et al 1995b, Williams et 
al 1996) around R = 22. The quasar and galaxy 
count estimates (converted to R band magnitudes 
using the colors of an f^, oc iy~^ power law) are 
plotted together in figure 4. We see that quasars 
are the more important foreground for R < 20. At 
fainter magnitudes, galaxies appear more impor- 
tant, but this may be misleading since they can 
often be identified using a third independent color 
and/or morphological information. 

Presently, most GRB error boxes have charac- 



6 



T ^ , , , I , , , ^ , , r 




R magnitude 

Fig. 4. — Cumulative number-magnitude relation 

for sources whose optical colors resemble GRB af- 
terglows. The solid lines shows quasar counts pre- 
dicted by models A and B of Malhotra & Turner 
(1995). The dotted line shows the predicted num- 
ber counts for galaxies with 0.4 < z < 0.6 us- 
ing the locally derived luminosity function from 
Loveday et al (1992). Counts of blue galaxies 
should be ~ 30% lower. The histogram shows the 
Canada-Prance Redshift Survey data for galaxies 
with 0.4 < z < 0.6 and with V — I colors in 
the afterglow range (Lilly et al 1995a; Le Fevre 
et al 1995; Hammer et al 1995). The total num- 
ber of galaxies around 22nd magnitude is a factor 
~ 10 above the galaxy counts shown. The plotted 
curves can be taken as an approximate guide to the 
number of potentially confusing sources expected 
for a given area searched and limiting magnitude 
attained. The exact number can depend substan- 
tially on additional factors, in particular photo- 
metric accuracy, spatial resolution, and the avail- 
ability of a third independent color to distinguish 
z ~ 0.5 galaxies from afterglows. 



teristic sizes of a few tens of square arcmimites. 
Thus, the brightest quasar expected in the error 
box has a magnitude R ~ 22, while the brightest 
galaxy having colors similar to an afterglow might 
be a magnitude brighter. Both this depth and 
the error box area are reasonably well matched to 
the capabilities of moderate sized (1 to 2.5 meter) 
ground-based telescopes. 

4.1. Optimizing the Observations 

To realize the promise of same-night afterglow 
identifications, a rapid and standard data analy- 
sis procedure is desirable, and could be aided by 
an observing procedure optimized for quick anal- 
ysis. Archival bias frames and flatfield frames for 
all filters could be useful, and precisely cospatial 
exposures in multiple filters (i.e., withoTit offsets 
and with telescope guiding) could ease measure- 
ment of object colors. 

An additional constraint on observing strategy 
arises because the target is likely to fade during the 
period of observations. If the multiple wavelengths 
cannot be observed simultaneously, it is then bet- 
ter to begin with the bluest filter, followed by the 
reddest, and only then move on to the middle fil- 
ters. This way, fading behavior will accentuate 
the expected color offset rather than reducing it. 
(The choice of bluest filter before reddest is driven 
by the possibility that a spectral break may move 
through the observed wavelength range during the 
course of observations.) 

The optimal choice of filter set will depend 
on local conditions. The relative difficulty of U 
band observations is approximately balanced by 
the large color-plane offsets obtained with U data 
(~ 0.7 mag from the stellar locus), which relax 
the signal to noise requirements considerably. 

4.2. Orphan Afterglow Searches 

If gamma ray bursts are highly coUimated, 
we expect to observe many more afterglows than 
GRBs (Rhoads 1997). The search for such "or- 
phan afterglows" is an important and largely 
model-independent way of testing GRB collima- 
tion. 

Color-based techniques will be a valuable ad- 
dition to orphan afterglow searches. Afterglows 
remain rare phenomena even in extreme coUima- 
tion scenarios, and in order to find one it will 
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be necessary to survey a large area of sky that 
will contain many quasars. Thus, orphan after- 
glow searches must combine variability with color 
tests to be reliable. When enough filters are ob- 
served simultaneously in a field where an archival 
first epoch is available, it will be possible to iden- 
tify new sources, test their colors for consistency 
with a power law spectrum, and start intensive 
followup of any viable afterglow candidates. This 
will allow much greater confidence in the identifi- 
cation of orphan afterglows than is presently pos- 
sible, since good spectra and light curves can then 
be used to distinguish real afterglows from vari- 
able stars, high redshift supernovae, active galac- 
tic nuclei, microlensing events, and other optically 
variable sources. 

As an illustrative example, consider Sloan Dig- 
ital Sky Survey (SDSS; see Knapp et al 1999). 
The SDSS southern strip will cover an area ~ 200 
degree^ to about 22nd magnitude in 5 broadband 
filters many (~ 45) times. The known GRB rate 
is ~ 3 per sky per day, and bright afterglows re- 
main above 22nd magnitude for a period ~ 4 days. 
Thus the southern strip has an effective coverage 
of 4 X lO^degree^ day. The northern SDSS survey 
will cover IC^ degree^ in a single imaging epoch. 
Fortunately, GRB afterglows would likely be tar- 
geted for spectra based on their quasar-like colors, 
yielding a second epoch. The northern survey thus 
provides another ~ 4 x lO^degree^ day. Both parts 
together would be expected to detect ~ 6 GRB af- 
terglows. If GRBs are tightly coUimated (opening 
angles < 10°), the expected number of orphan af- 
terglows rises to > 25. By using the color space 
properties to identify these afterglows, such a sur- 
vey can place strong limits on the collimation an- 
gle and hence energy requirements of gamma ray 
bursts. 

5. Conclusions 

Color-based searches will allow faster identifi- 
cation of many gamma ray burst afterglows. In- 
tensive followup of these afterglows can then begin 
earlier and at brighter flux levels, yielding a larger 
sample of bursts for which spectroscopic redshifts 
and well measured spectral energy distributions 
(including roughly synoptic data at X-ray, opti- 
cal, near-infrared, subinillimetcr, and radio wave- 
lengths) can be obtained. Such data sets are cru- 



cial to detailed studies of afterglow physics. Ad- 
ditionally, single-epoch color-based searches may 
yield a more easily characterized statistical sam- 
ple of afterglow detections and nondctections than 
is currently available. Such samples will allow reli- 
able studies of the population properties of gamma 
ray burst afterglows, and thereby open the way for 
new insights into the nature of gamma ray bursts. 
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